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Abstract 

The original Calogero and Sutherland models describe N quantum particles 
on the line interacting pairwise through an inverse square and an inverse sinus- 
square potential. They are well known to be integrable and solvable. Here we 
extend the Calogero and Sutherland Hamiltonians by means of new interactions 
which are PT-symmetric but not self adjoint. Some of these new interactions 
lead to integrable PT-symmetric Hamiltonians; the algebraic properties further 
reveal that they are solvable as well. We also consider PT-symmetric interac- 
tions which lead to a new quasi-exactly solvable deformation of the Calogero 
and Sutherland Hamiltonians. 



*yves. brihaye@umh.ac.be 
^ nininaha@yahoo . fr 



1 Introduction 

The Calogero and Sutherland models [1, 2] (CSM in the following) and various of their 
extensions have received a considerable new impetus of interest in the last years. We 
mention for example the recent preprint [3] where CSM are applied in the framework 
of conformal field theory. From the beginning, the CSM are know to be integrable 
both on the classicall and quantum levels. In the framework of quantum systems, 
the distinction between integrable and exactly solvable models can further be made, 
according to the lines of Ref. [4] where it was shown that quantum CSM possess both 
properties: integrable and exactly solvable. Off course the notion of exactly solvable 
operators constitutes a particular case of Quasi-Exactly-Solvable (QES) operators [5] 
for which only a finite part of the spectrum can be computed algebraically. The 
possibihty of having QES extensions of the Calogero and Sutherland Hamiltonians is 
a natural question which was adressed e.g. in [6, 7]. 

Recently again, it was recognized that the self-adjoint property is not necessary 
for an operator A to have a real spectrum. Instead the weaker condition of invariance 
of A under the combination of parity P and time-reversal T symmetries leads to a 
spectrum which is either real or composed of pairs of complex conjugate numbers [8]. 
It is therefore natural to consider non-self adjoint but PT-symmetric Hamiltonians 
describing N-particles quantum systems and to construct operators which are (i) 
integrable or (ii) exactly solvable or (iii) quasi-cxactly solvable. Examples of such 
operators were studied in [9] and reconsidered more recently [10] 

Off course the problem of classifying all PT-symmetric (quasi)- exactly integrable 
operators is vast, one possible way to approach it [9, 10] is to extend the well estab- 
lished CSM by suitably chosen non self-adjoint but PT-symmetric extra terms and 
to study which algebraic properties of the CSM are preserved or spoiled by the new 
terms. In this paper we consider more general types of PT-symmetric terms restricted 
by the following properties (i) translation invariant, (ii) first order in the derivatives, 
(iii) fully symmetric under the permutations of the N degrees of freedom. In the three 
cases we put the emphasis on the (quasi) solvability of the deformed model. We will 
exhibit two extensions which preserve the exact solvability and another one leading 
to a new type of QES N-body Hamiltonian. 

2 PT-invariant Hamiltonians: rational case 

We consider the quantum N-body Hamiltonians of the form 



H = Heal + HpT 



(1) 



where Hcai denotes the rational Calogero[l] Hamiltonian : 
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6Hq + jHi 
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f2 = ^{xj - Xkf , ai = J2xj , (5) 

jj^k j=l 

For later convenience, we split the traditional harmonic part into the translation- 
invariant piece (f2) and the "center of mass coordinate" (7i; 

Setting e = the above Hamiltonian is purely translation invariant, the standard case 
is recovered for e = 1. The part HpT is not self-adjoint but symmetric under the PT 
reflexion i ^ —i , Xk ^ —Xk , Pk ^ Pk ipk = f ^f^)) it extends the interaction Hq 
considered recently in [10]. In order to reveal the solvability property of H we first 
perform the standard (often called "gauge" ) transformation 

H^-^-^H-^gr , '^gr = (3''E (7) 

with 

(3=Yl{xj-Xk) , E = exp(-|f2 - %D (8) 

j<k 

After some algebra, the new Hamiltonian H is obtained. The parameters u, a, a' 
entering in "ifgr are fixed according to 

a'N 

g^y'^ 25) , uj"^ = a'^N + 2a^N , uj = -— (9) 

in such a way that the terms with the highest power in Xi as well as the most singular 

terms in Xi disappear in H. The condition on g, p, S is jiist identical as in the case of 
Ref. [9] and lead to the same restrictions on these constants for the eigenvectors to be 
nonsingular : (i) (5 > — | , > 51 > —{S + and 51 > for arbitrary values of 5. 



2.1 Exact solvability of H 

The next step consists in writing H in terms of variables which encodes the full 

symmetry of the problem under permutations of the N degrees of freedom. We use 
here the variables introduced in [4] which we recall for completeness : 

(Ji(x) = Xi + X2 + . . . + xn- (10) 

a2{x) = X1X2 + xiXs + . . . + xn-iXn (11) 

(^si^) = E ^i^j^k , ■■■ (12) 
i<j<k 

crwix) ^ X1X2X3. . .xn (13) 
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and ior k — 2,3, . . . , N we further define tlie translation invariant variables 

Tkix) = ak{y) Vj = Xj - -^cTi (14) 

In particular, we have f2 = —4Nt2. The Laplacian operator once expressed in terms 
of these new variables reads [4] 

N Q2 

3=1 ^-^3 

Q2 N Q2 N 

= N^,+ T.A>.^ + Y.B,-I^, (16) 

^^1 i,fc=2 OTjOTk ^^3 

l>max{l,k—j) 

{N-j+2){N-j + l) 
-^3 j;^ ^3-2 (l^j 

where the following conventions have to be understood 

To = 1 Ti = T_p = tn+p = for p > (19) 

The PT invariant interactions Hq, Hi defined in Eq.(3) can also be expressed in terms 
of the new variables by means of the following identities : 

N Q 

D =1:^JqZ- (20) 

j=l ^Xj 

d ^ d 

= ^i^ + E^^^^ (21) 

5(71 

.a.A + l, (22) 

I N o 

^0 ^--Y.iN-j + 2){N-j + l)T,.,— (23) 

j=2 3 

Hi ^ND (24) 
The final form of the operator H reads then 

H =^A + {u-5)Ho 



+u;^l + ^-^D + u;eai— + Eo (25) 
£^0 =«^^ ^(i/7V + i) + _- + i.^^- L-Sa^^ ^ (26) 
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The solvability of the operator H (and therefore of H in (1)) is revealed by the 
observation [4] that this operator (25) preserves the infinite flag of vector space 

N 

Vn = span«V2"^ ...T^\Y.^k<n), ne-H (27) 

k=l 

The eigenvectors of H can therefore be constructed by considering the restriction of 
^ to V„. As far as the construction of the eigenvalues is concerned, it can be realized 
that the restriction of H to Vn leads effectively to a lower triangular matrix, the 

diagonal elements of which are issued from the matrix elements corresponding to the 
second line in Eq.(25). Off course the triangularity of this matrix is strongly related 
to the order adopted to enumerate the basis elements r2 ^ . . . r^^ . The relevant order 
can be set as follow 

< rijsi < n , 

< nAr_i <n — un, ... , 

< 712 < n — nN — nN-i ■ ■ ■ — (28) 

The eigenvalues can therefore be determined in terms of the diagonal matrix 
elements by using the monomials defined in Eq.(27) as a basis. The spectrum reads 

Ki,n2,...,n;v = ^0 + ujem +UjJl + ^^(2^2 + Srig + . . . Nun) (29) 

The spectrum of the Calogero model is off course recovered with all its degeneracies 
in the limit 7 = 0; the term Hq does not affect these degeneracies. However for 7 7^ 
several degeneracies, but not all, are hfted by the Hi interaction. The eigenstates 
which remain degenerate are those with J2k=2 ^''^k — J2k=2 ^'^fe- 

2.2 Quasi-exactly Solvable case 

In this section we consider the Hamiltonian 

H^H^i + eY, (^i - xjf^ + eJ2 - ^j)' (so) 

where H^ai is supplemented by a PT-invariant term (with coupling constant 6) as 
well as a new term of the potential. This term is anharmonic, being of degree 4 in 
the position, and characterized by a coupling constant 9. Performing the change of 
function by means of ipgr leads to 

H ^ -^A + uHq + cuD 

,~ 9a. ,4 ,a'^N uP' ^ ^-^ f^^ /o-,x 

+ (^ (^^ - ^i) + 4iV^2(^ - 47V " " ^^^^ 
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If the new coupling constant 6, 9 are chosen in such a way that 

then the operator H preserves a finite dimensional vector space which we call V^. 
As so, it is quasi-exactly solvable [5]. The vector space Vm is characterized by the 
condition 

^ d ~ ~ 

T2{m - D)Vm ^ T2{m - ^ kTk-^)Vm <^ Vm (33) 

k=2 ^'^k 

A little algebra shows that m can not be arbitrary. As an example, if = 3 the 
vector space V^n has the form 

Vm = VC^) ® r|P(| - 3) © rtvC^ - 6) ® . . . (34) 

if m is even (m > 2) and 

Vm = r3P(^— ) © t|^(^— ) ® t|^(^— ) ® ■ ■ ■ (35) 

if m is odd (m > 5). The sum off course runs with k as long as (m — 6k) /2 (resp. 
(m — 3 — 6/c)/2) is positive for m even (resp. odd). Remark that we have the following 
inclusions 

Vm C Vm/2 , Vm C V(m-3)/2 (36) 

and that the vector spaces Vm are preserved separately by the operators A, Hq, D. In 
order to illustrate these results, we constructed the eigenvalues of the Hamiltonian in 
the case A^ = 3, m = 6, for which the vector space Vq has five dimensions.The couphng 
constant 9 of the quartic interaction in (30) is then adjusted in terms of 9, ui, v and 
m by Eq (32). Note that Changing the value of m results in another value for a and 
therefore another value for 9. The eigenvalues E cannot be computed explicitely in 
terms of the coupling constants v,uj,9 but, considering the PT-symmetric term as a 
perturbation of the Calogero Hamiltonian leads to the following result 

£;o,o = + 27(1/ -l)- + 0(0 

Q 

^1,0 =2a; + 9(z/-5)- + 0(^2) 
E2,o ^ Aco - 9{u + 1)- + 0{9^) 

UJ 

Es^o = 6a;-27(i/-3)- + 0(^2) 

Eo,2 =6u; (37) 

Note that the only degenerate energy levels which occur in the diagonalisation of the 
restricted to V3 corresponding to m2 = 3, ms = and m2 = 0, ms = 2. These vectors 
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are included in the subspace Vq and, again, this new term hft the degeneracy. The 
six extra vectors corresponding to V3 are not accessed algebraically in the presence of 
the new term. We also solve numerically the eigenvalue equation for generic values 
of u, 9 and found that the spectrum remain real. 



3 PT-invariant Hamiltonians: trigonometric case 

Here, we consider an extension of the Sutherland Hamiltonian Hgu in the form 

H ^H,u + HpT + V (38) 

where the extra PT symmetric part is chosen of the form 

HpT ^5Ho + -fHi (40) 
^5^cot(^^^)( ^ 



i<j 2 dxi dxj 

+7Esin(xi-^i)(T|r-^) (41) 

which provides a natural generalisation of Eq.(3) to the trigonometric case. The 
Sutherland inverse sine-square potential is also supplemented by the following terms 

V =e,j: cot(:^)cot(^i^) 

+^3 E cos(^^)cos(^i±^i^) (42) 

where the ^2-term involves two body-interactions (like the original Sutherland poten- 
tial) while the Q\ and ^3-terms involve three body-interactions. Similarly the previous 
case, it is convenient to "gauge rotate" of B. in (38) by using the ground state ■^pr of 
the Sutherland model i.e 

v^.. = (n«M^^)r ^^+^=^ (43) 

3<k ^ 
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In this purpose, the following identities are usefuU 

o-lrro rr , y 1 N{N - l) 

+ E cot(^^)cot(^^) (44) 
P-'H,P ^H, + Y^cos'{^^)+ Y: cos(^i^)cos(^i±^i^)(45) 

i<j i¥=j¥=k 

The gauge rotation of the full Hamiltonian (38) by ■0sr leads to the following 
equivalent operator 

h ^ -2f3-''Hf3'' (46) 
= A + (i/ - 2S)Ho + {-2j)Hi + Eo (47) 

where i?o = i''^N{N'^ — 1)/12 denotes the energy of the ground state. In order to 
obtain the above formula, the various coupling constants are choosen according to 

^ = i.(i.-l) + 5 , , ^2 = ^3^-7 (48) 

in such a way that the singular terms and the non derivative parts of the operator h 
cancel. 

3.1 Exact solvability 

The next step to reveal the solvabihty of the Hamiltonian (46) is to use the following 
change of variables 

= e*^i + e^^2 + . . . + e'^^ (49) 
jyi = e'v^ + e'y' + ... + e'y^ (50) 

_ g«(2/i+J/2) _|_ gifei+'ya) gifeiV-i+J/Jv) ^5]_-j 



V2 



1 ^ 



77JV-1 =771 > Vi^ Xi- — ^Xk (52) 

fe=i 

The different operators can then be expressed in terms of the new variables namely 
^ifc = W + E {k-j- 2l)rij+irik-i (54) 

i>maa;(l,fe— 

i^o =-E^(^-0^^^ (55) 
1=1 ^'ii 



8 



Considering first the Hq interaction only (7 = ^2 = ^3 = 0), we observe that this 
PT-symmetric term does not spoil the solvability of the Sutherland Hamiltonian [4] 
because both A, Hq preserve the infinite flag of vector spaces Vn{v) (see (27) for the 
definition) as easily seen from Eq. (46). It is, however, worth to stress that the term 
Hq leads to solvability only if the extra three-body potential corresponding to the 61- 
term in Eq. (42) is supplemented to the original two-body Sutherland potential. This 
contrasts with the rational case where the addition of the term Hq to the Calogero 
Hamiltonian leads to another solvable model without any extra pieces in the potential. 
The interaction Hi is investigated in the next section. 



3.2 Quasi Exact solvability 

We could not evaluate Hi in terms of 77 for generic cases of N but, for N=2,3, we find 
respectively 

Hi =^(^1-4)^1^ ' 77i = 2cos^^ (56) 

Hi =^(r/iV-2%'-377i)^ + i(?72V-2r^?-3r/2)^ (57) 

rji ^ e'^' + e'y' + e'^' , r]2 ^ v*i (58) 

We see that, for 7 7^ 0, the Hamiltonian h is hkely nor solvable neither quasi-exactly 
solvable as suggested by the form of Hi for N=3, (see Eq(57)). This operator obviously 
does not preserve any of the Vn{i]) (see (27)). 

For N=2, the operator (38) is nevertheless QES provided 62 = — 27n, n G IN in 
(42). Indeed, the "gauge rotated" operator reads in this case 

which can be easily checked to preserve the vector space 

Vn = span(r/r , Vi''' , Vi'' ■ ■ ■) C V„ (60) 

Accordingly, (n-|-2)/2 (resp. {n + l)/2) eigenvectors can be constructed algebraically 
for n even (resp. odd). 

Once more we notice that, in spite of the fact that (40) constitutes the natural 
counterpart of Eq. (3) for the trigonometric case, the solvability is preserved only for 
the interaction Hq; in contrast to the rational case, the Hamiltonian Hi spoils the 
solvability of the model for N > 2. 



4 Conclusion 

We have considered several extensions of the Calogero and Sutherland models by 
means of PT-symmetric terms which are translation invariant, completely symmetric 
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and involve first derivatives. The algebraic structure of these extended Hamiltonians 
is revealed by (i) the change of function involving a common ground state function 
ipgr, (ii) the change to appropriate symmetric coordinates and(iii) the construction 
of some vector spaces of polynomials preserved by these operators. The completely 
solvable character of the Calogero Hamiltonians is preserved for the terms called Hq 
and Hi. Another choice of interaction, involving a PT-symmetric term supplemented 
by an extra anharmonic potential, leads to a new Quasi-Exactly-Solvable extension 
of the Calogero model. Note that this QES extension differs from the one obtained 
first in [6] where a potential depending of the variable T2 was added. 
Applying the same kind of methods to the trigonometric (or Sutherland) model re- 
veals that some PT-invariant terms lead to completely solvable extended models only 
provided a suitable extra three-body interaction is added to the potential. For an- 
other natural PT-symmetric term, the solvability is spoiled even at the price of adding 
new terms to the Sutherland original potential. Only for = 2 the corresponding 
operator turns out to be quasi-exactly solvable. 
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